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1. Introduction 

The frequently discussed model of twelve fermions in the fundamental representation of the 
SU(3) color gauge group remains controversial with continuing recent efforts from several lattice 
investigations [1-13] where more extended references can be found for the earlier history of the 
model. The focus of these lattice studies is to establish conformality or chiral symmetry breaking 
in the bulk phase close to the continuum limit. At finite cutoff, two different strategies can be used 
to deal with finite volume dependence while probing the conformal and #SB hypotheses of the 
bulk phase. The first strategy extrapolates the spectrum to infinite volume at fixed fermion mass 
m where the leading finite size correction is exponentially small and assumed to be determined by 
the lowest mass with pion quantum numbers. From the mass spectrum of the infinite volume ex- 
trapolation the mass deformed conformal scaling behavior can be probed and compared with ^SB 
behavior when the fermion mass is varied in the infinite volume limit. This limit requires very large 
volumes and it is difficult to reach with existing computational resources. Probing the conformal 
hypothesis without intrinsic scale, the second strategy takes full advantage of the conformal FSS 
behavior when approaching the m = critical surface at fixed finite size L. Different from the first 
strategy, the finite volume corrections are not exponentially small and a much larger data set can 
be analyzed closer to the critical surface. Results from the conformal FSS analysis are presented 
here, significantly extending what was reported earlier [1]. 

2. The algorithm and the data base from the simulations 

We have used the tree-level Symanzik-improved gauge action with /3 = 2.2 gauge coupling 
and with two levels of stout smearing in the staggered fermion action [14]. The standard HMC 
algorithm was used with multiple time scales and the Omelyan integrator. Our error analysis of 
hadron masses used correlated fitting with double jackknife procedure on the covariance matrices. 
The time histories of the fermion condensate, the plaquette, and correlators are used to monitor 
autocorrelation times in the simulations. 

Extending earlier work [1], we have new simulation results at jS = 2.2 in the fermion mass 
range m = 0.002 - 0.025 at lattice volumes 20 3 x 40, 24 3 x 48, 28 3 x 56, 32 3 x 64, 40 3 x 80, and 
48 3 x 96. The extended data base now has the m = 0.002 — 0.035 fermion mass range with a span 
close to a factor of twenty. The new low fermion mass range m = 0.002, 0.004, 0.006, 0.008 is used 
in the conformal FSS analysis which over the full set extends the lattice pion correlation length 
from 2.5 to 20 in the infinite volume limit. Results from the two lowest masses at m = 0.002, 0.004 
are not included in the current analysis and will be reported later. For further control on finite 
volume dependence, large 48 3 x 96 runs were continued to two thousand trajectories at m = 0.01 
and m = 0.015. Four runs were further added at 40 3 x 80 with m = 0.01,0.15,0.02,0.025. The 
new and updated data sets were subjected to conformal FSS analysis and #SB tests of the (xjfY) 
chiral condensate. Preliminary analysis of the results already appeared a few months ago [2]. 

3. The bulk phase diagram 

The bulk phase structure of the model remains controversial, particularly the critically im- 
portant weak coupling phase. In addition to our spectroscopy and conformal FSS runs, we ran 
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extensive scans at various fixed volumes and fixed fermion masses to explore the bulk phase struc- 
ture. The bare coupling /3 was varied over a large range starting from very small /3 values deep in 
the strong coupling regime to the weak coupling phase at /3 = 2.2 where the conformal and #SB 
analyses were done. Fermion masses m = 0.007,0.01,0.02 were used in the scans with spatial 
lattice sizes L = 8, 12, 16,20,24,32 running a large densely spaced set in the important and much 
discussed intermediate region in transit from strong coupling to weak coupling. These scans were 
also extended to Nf = 2,4, 6, 8, 10, 12, 14, 16 flavors. Two representative scans of the bulk behavior 
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Figure 1: On the left, scans of the phase diagram by monitoring the chiral condensate are plotted as a function of /3 at 
two different fermion masses for Nf = 12. The schematic bulk phase diagram is sketched on the right, a plot suggested 
from findings in [9] . 

of the chiral condensate {yy} are shown in Figure 1 as we vary j3 from strong to weak coupling. 
Three distinct regions emerge at fixed volume and fixed fermion mass showing strong coupling 
behavior for j3 < 1.4 with a large chiral condensate, an intermediate phase for 1.4 < j8 < 1.8 with 
a sudden drop in {yy}, and a weak coupling phase for /3 > 1.8 with a further drop in (yy). We 
also observed another distinct feature of the intermediate phase (first seen in [9]) where the pion 
correlator exhibits oscillations associated with the appearance of a non-degenerate parity partner 
state. This structure does not exist in the weak coupling phase. A similar pattern of phases was also 
seen in scans at Nf = 8. Our physics simulations and conformal FSS hypothesis tests were done 
well inside the weak coupling phase at /3 = 2.2 without influence from the intermediate phase, as 
indicated in Figure 1. 

A similar phase structure has been observed independently by three different groups [4,9, 12]. 
The newfound order parameter of broken shift symmetry in the intermediate phase is an interesting 
additional development in the study of the esoteric intermediate phase when HYP smearing is used 
in the staggered fermion action [9]. The intermediate phase only exists in a finite interval of the 
lattice gauge coupling for small enough fermion masses, as schematically sketched in Figure 1. 
The real interest is, of course, in the nature of the weak coupling phase. Based on axial U(l) 
symmetry considerations, arguments were presented in [12] in favor of conformal symmetry in 
the weak coupling phase. This argument was criticized in [9] based on new details of the broken 
shift symmetry with chiral symmetry restoration they discovered at zero temperature in the bulk 
intermediate phase. 

Weak coupling results on the Polyakov loop, the chiral condensate, and spectroscopy were also 
presented in [9] as indications of conformal symmetry in the weak coupling phase. The blocked 
Polyakov loop was reported to jump from zero to a large finite value in crossing to the weak 
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coupling phase [9] . A confining potential was reported in the intermediate phase with broken shift 
symmetry which turned into a Coulomb potential without a string tension in the weak coupling 
phase [9]. It was also asserted that the observed chiral condensate and the related Dirac spectrum 
show the recovery of exact chiral symmetry in the massless fermion limit of the weak coupling 
phase consistent with observed degeneracies of parity partners even at finite fermion masses. In 
conclusion, conformal behavior was suggested for the bulk weak coupling phase. 

The results reported in [9] suggesting a chirally symmetric deconfined conformal phase are 
in disagreement with what we found earlier [1] and in the extended new analysis [2,3]. Using 
lattice volumes several times larger than the simulations in [9] we find a vanishing Polyakov loop 
at zero temperature in the weak coupling phase and a confining potential at a pion correlation length 
which is significantly larger and consequently more relevant than the related findings in [9] . We also 
observe the splitting of parity partners at finite fermion mass. Consequently, and differing from [9], 
our findings in large volumes are consistent with a chirally broken weak coupling phase. We turn 
now to conformal FSS tests of the conformal hypothesis in the weak coupling phase gathering 
further evidence towards more definitive conclusions. 

4. Conformal finite size scaling analysis 

The expected leading FSS form for any mass M in the spectrum, or for the decay constant F n , 
scaled with the linear size L of the spatial volume, is given by a scaling function L-M = f(x) where 
x = L- m 1 / 1 ^ is the conformal scaling variable. The scaling form sets in close to the critical surface 
for small m values. The scaling functions f(x) can depend on the quantum numbers of the states but 
the scaling variable is expected to have the same form with identical y exponent in each quantum 
number channel [15-19]. In sub-leading order there are conformal FSS scaling violation effects 
which are exhibited as a combined cutoff and L-dependent leading correction with the modified 
form L-M = f(x) + L~ ( °g(x) where the scaling correction exponent CO is determined at the infrared 
fixed point (IRFP) g* of the /3 -function as CO = j3'(g*). This assumes that the mass deformation 
away from the critical surface is the only relevant perturbation around the IRFP. The leading scaling 
correction term close enough to the critical surface dominates any other corrections which are 
further supressed by additional inverse powers of L. To detect the leading scaling violation effect 
requires high precision data with fits to scaling functions f(x) and g(x) and the critical exponent 
CO. 

4.1 Conformal FSS fitting procedure with restricted scaling functions 

We applied conformal FSS theory to our data sets in the fermion mass range m = 0.006 — 0.035 
with lattice sizes ranging in the fits from 20 3 x 40 to 48 3 x 96. Two different FSS fitting procedures 
were applied. In the first procedure, we defined a scaling function f(x) for each mass M with five 
independent fitting parameters. The fitting function f(x) was divided into two regions separated at 
the joint x = x cut . Different forms were chosen on the two sides of x cut for the expected conformal 
behavior. For large x > x cut , the function f(x) = c\x + c exp (cix)~ 1 ' 2 exp(— c\x) with parameters 
c\ and c exp describes the L-independent limit M ~ cim 1//1+r at fixed m and L — > oo. The c exp 
amplitude sets the size of the leading small exponential correction from the finite volume effect of 
the lightest Goldstone pion state wrapping around the spatial volume. Since /(0) = cq is expected 
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conformal FSS (Goldstone pion - PCAC channel) 
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Figure 2: Conformal FSS fits in four different quantum number channels. The fits are performed in each 
channel separately. Since the y values vary considerably from channel to channel, a simultaneous global fit 
to the combined channels with the same / exponent, as required by conformal FSS theory, is bound to fail. 



from conformal FSS with some power corrections at small x, we applied the simple ansatz f{x) = 
Co + c a x a for x < x cut (a more general polynomial function in the small x region will not change the 
conclusions from the fits). From the fit to the PCAC Goldstone pion channel the parameter c % = ci 
was determined and used as input in the exponential terms of the other channels with exp(— c K L). 
The critical exponent y was included among the five fitting parameters, in addition to cq, c\, c exp , 
nn&Xcut. 

The composite particle masses in several quantum number channels can be reasonably fitted 
with conformal scaling functions f{x) as shown in Figure 2 but the values of the critical exponent 
7 are incompatible across different channels. The required global conformal FSS fit will fail with 
a single exponent y across all quantum numbers. In the fits for F n in the PCAC pion channel 
we only kept four parameters because the asymptotic form with exponentially small correction 
was zero within error. Actually, the data of F K did not allow a successful conformal fit with any 
shape chosen for its scaling function f{x) which looks very different from the scaling functions of 
composite particle masses. The unexpectedly curious behavior of the F n data set against conformal 
FSS remains unresolved. Our dataset is not accurate enough to successfully resolve subleading 
corrections in the fits. 
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4.2 Generalized FSS fitting procedure with spline based general B-form 

Following a new fitting strategy, we investigated if the failed global conformal FSS analysis 
can be attributed to restrictions on the conformal scaling functions f(x). The restrictions were 
manifest in the physics-motivated fitting procedure we applied above. Our new general approach 
is different from [20,21] but addresses related issues. We developed a general least-squares fitting 
procedure to the scaling functions using the B-form of spline functions [22] without additional 
restrictions. In this procedure, the function f(x) is described by piece-wise polynomial forms con- 
structed from spline base functions with general coefficients in overlapping intervals of the scaling 
variable x which depends on the exponent y. The shape of the B-form can be changed without 
limitations by increasing the number of base functions and the number of scaling intervals in x 
bracketing the overlapping data range. The details of this new analysis will be reported elsewhere. 



conformal FSS fit with spline based B-form (Goldstone pion) conformal FSS fit with spline based B-form (cRho2 channel) conformal FSS fit with spline based B-form (F in PCAC channel) 




x= Lmy m x= Lmv m ' x= Lmv m ' 



Figure 3: Conformal FSS fits using spline based B-forms in three different channels. The fits are performed in each 
channel separately with the question mark on y indicating difficulties of error estimates in bad fits of F K . 

Our fitting procedure in its setup requires two steps. In the first step, for any fixed choice 
of the exponent y, the best fitted function f(x) is determined in spline function B-form from the 
minimization of the weighted % 2 expression. According to a general algorithm, the x-range of the 
data set is divided into intervals separated by internal knots and adding end point knots for B-form 
spline construction. The number of coefficients is determined by the number of knots and the order 
of the spline polynomials of the sub-intervals. The weighted % 2 sum is minimized with respect 
to the coefficients of the base functions in the B-form. This will produce the best fit for fixed y 
with a minimized % 2 sum which will depend on y. In the second step, we minimize the % 2 sum 
with respect to y to determine the best fit of the critical exponent. The one-a confidence interval is 
determined from the variation of the % 2 sum as a function of y. 

In Figure 3 we show three typical fits for illustration. The fit to the Goldstone pion in the 
PCAC channel improved as expected, with considerable increase in the error. The tension across 
channels decreased, as illustrated by comparison with the rho-channel fit, but the fit to F n remained 
unacceptable. With the extended data set we are unable to reproduce results in [20,21] which used 
tables from our earlier limited subset of data [1] to argue in favor of consistency with the conformal 
phase. It is important to emphasize that we have not reached definitive conclusions about the failure 
of conformal tests. As we stated earlier [1], we have not analyzed yet the leading scaling violation 
effects and did not investigate if the good scaling form in separate quantum number channels can be 
explained in the chirally broken phase by strongly sqeezed wave function effects. In disagreement 
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with [21], conformal FSS based analysis of the spectrum and related sum rules on moments of the 
correlators we have been developing are deep renormalization group based probes of the conformal 
phase. As explained in our forthcoming publication, we remain skeptical about the fitting procedure 
followed in [20] with efforts to rescue the conformal interpretation. The issues are not settled and 
ultimately will be decided in more definitive analyses. 
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